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Abstract 

The Hofstadter butterfly is a quantum fractal with a highly complex nested set of gaps, 
where each gap represents a quantum Hall state whose quantized conductivity is char¬ 
acterized by topological invariants known as the Chern numbers. Here we obtain sim¬ 
ple rules to determine the Chern numbers at all scales in the butterfly fractal and lay 
out a very detailed topological map of the butterfly. Our study reveals the existence 
of a set of critical points, each corresponding to a macroscopic annihilation of orderly 
patterns of both the positive and the negative Cherns that appears as a fine structure in 
the butterfly. Such topological collapses are identified with the Van Hove singularities 
that exists at every band center in the butterfly landscape. We thus associate a topo¬ 
logical character to the Van Hove anomalies. Finally, we show that this fine structure 
is amplified under perturbation, inducing quantum phase transitions to higher Chern 
states in the system. 


Discovered by Belgian physicist Leon Van Hove in 1953, Van Hove singularities 
are singularities in the density of states of a crystalline solidQ. These singularities are 
known to be responsible for various anomalies provided Fermi level lies close to such a 
singularity. Electronic instabilities at the crossing of the Fermi energy with a Van Hove 
singularity in the density of states (DOS) often lead to new phases of matter such as 
superconductivity, magnetism, or density waves El. 
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A two-dimensional electron gas (2DEG) in a square lattice provides a simple exam¬ 
ple of Van Hove singularities in the energy dispersion of a crystal. For a tight binding 
model of a square lattice , the energy dispersion is given by, 

E = — 2J[cos k x a + cos k y a] (1) 

Here k = (k x ,k y ) is the wave vector in the first Brillouin zone and a is the lattice 

spacing of the square lattice and J is the nearest-neighbor hopping parameter which 
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defines the effective mass m e of the electron on the lattice by the relation J = 2y ^ q2 . 
This single band Hamiltonian has band edges at E = ±4 J. It can be shown that the 
density of states (DOS) at the band edges approaches a constant equal to • How¬ 

ever, it diverges at the band center as DOS « lj ^~. Such a divergence is an example of a 
Van Hove singularity. Figure [T] shows the energy contours in (k x ,k y ) plane, where the 
almost free-electron concentric circles are transformed into a diamond shape structure 
that corresponds to saddle points in the energy surface. We note that the lattice struc¬ 
ture is essential for the existence of Van Hove singularities. Van Hove singularities 
have been given a topological interpretation in terms of a switching of electron orbits 
from electron like to hole like 01. 

In this paper we investigate the Van Hove anomalies of a 2DEG in transverse mag¬ 
netic fields. Such system describes all phases of non-interacting electrons as one varies 
the chemical potential and magnetic field. The phase diagram, known as the Hofstadter 
butterfly (TOl represents various quantum Hall states, each characterized by a quantum 
number, the Chern number, that has its roots in the nontrivial topology of the underly¬ 
ing Hilbert space Q. The key result of this paper is the topological characterization of 
Van Hove singularities that are nested in the hierarchical pattern of the butterfly spec¬ 
trum. We show that in the two-dimensional energy-flux space, every vicinity of a Van 
Hove consists of interlacing sequences of positive and negative Chern numbers that 
collide and annihilate at the Van Hove singularities. In other words, Van Hove singu¬ 
larities induce a topological collapse in the quantum fractal of the Hofstadter butterfly. 
We calculate Chern numbers in the neighborhood of Van Hove singularities, facilitated 
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Figure 1: (a) Contour plot of the energy E in the (k x ,k y ) plane, illustrating the saddle character of the band 
center for a 2DEG on a square lattice, (b) Shows the corresponding band as a function of k x . (c) shows the 
spectrum for small magnetic flux </> = 0.1. Landau levels correspond to the horizonal flat bands. In (d) we 
show a blow up of c) near the band center, illustrating the deviation from the Landau level picture near the 
band center that hosts a Van Hove singularity. In (b),(c) and (d), the colours represent different values of k y . 


by simple rules that we derive for determining the entire topological map of the but¬ 
terfly fractal at all scales. Our analysis begins with a geometrical approach, based on 
simple number theory, that sets the stage for determining the Chern numbers of all the 
gaps and its associated fine structure. It is the orderly patterns of topological integers 
that characterize the fine structure that gets linked to the Van Hove anomalies of a two 
dimensional crystalline lattice in a magnetic field. In other words, we correlate the 
macroscopic nature of the topological collapse from which emerges a reincarnation of 
Van Hove in a quantum fractal made up of integers characterizing the quantum Hall 
conductivity. 


A very recent study of the 2DEG when subjected to a weak magnetic field®, re¬ 
vealed the importance of Van Hove singularities in inducing changes in localization 
characteristics of the system. In a continuum system, that is, in the absence of any 
lattice structure, the magnetic field B introduces a magnetic length l b = yj 2 Eb » re i n_ 
carnation of the cyclotron radius of the corresponding classical problem. In this limit, 
the energy spectrum consists of equally spaced harmonic oscillator levels known as 
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the Landau-levels. Interestingly, in a lattice with weak magnetic flux, the Landau level 
picture breaks down near the band center as illustrated in the Fig. 0. The insight into 
this clustering of the band fragments near the center can be gained by expanding the 
energy E near k x -\-k y = ± 7 r and k x — k y = ± 7 r, where the energy depends linearly on 
the wave vector. This is in sharp contrast to the quadratic dependence of energy near 
the band edges that leads to simple harmonic levels, namely the Landau levels. 

The model system that we study here consists of (spinless) fermions in a square 
lattice. Each site is labeled by a vector r = nx + my, where n, m are integers, x ( y ) 
is the unit vector in the x ( y ) direction, and a is the lattice spacing. The tight binding 
Hamiltonian has the form 

H = - J*E l r + *>< r l - Jy E l r + vy 2m HA + h.c. (2) 

r r 

Here, |r) is the Wannier state localized at site r. J x ( J y ) is the nearest neighbor hop¬ 
ping along the x ( y ) direction. With a uniform magnetic field B along the 2 direction, 
the flux per plaquette, in units of the flux quantum &o, is (j> = —Ba 2 /Q o- Field B 
gives rise to the Peierls phase factor e z27m< ^ in the hopping. 

Within the Landau gauge, the above Hamiltonian has been engineered in cold atom 
experimentsQ. In this case, the vector potential is given by A x = 0 and A y = —<\>x 
resulting in a Hamiltonian that is cyclic in y. Therefore, the eigenstates of the system 
can be written as 4/ n?m = e' lkyl7l 'ip n where ^ n satisfies the Harper equation (TOl 

e^Vn+i + e-^Vn-i + 2A cos(2tt n</> + k y )^ = Eijj r n . (3) 

Here n (m) is the site index along the x ( y ) direction, A = J y /J x and Vn+q = Vn, 
r = 1, 2, ...q are linearly independent solutions. In this gauge, the magnetic Brillouin 
zone is —7r jqa < k x < n/qa and — 7 r < k y < tt. 

At the rational flux 0 = p/q, where p and q are relatively prime integers, the en¬ 
ergy spectrum has q — 1 gaps. These spectral gaps are labeled by two quantum numbers 
which we denote as a and r. The integer a is the Chern number , the quantum number 
associated with Hall conductivity (71 and r is an integer. For a Fermi level inside each 
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energy gap, the system is in an integer quantum Hall state(6) characterized by its Chern 
number a which gives transverse conductivity 171 C xy = cr^-. 

The quantum numbers (cr, r) satisfy a Diophantine equation (DEi fl4l . that applies 
to all 2DEG systems that exhibits magnetic translational symmetry, 

p = (\)G + r (4) 

where p is the particle density when Fermi level is in the gap. For a given p and 0, 
there are infinity of possible solutions for where [cr, r] are integers, given by, 

[cr, r] = [cr 0 - nq , r 0 + np\ (5) 

Here cr 0 , tq are any two integers that satisfy the Eq. 0 and n is an integer. The 
quantum numbers cr that determine the quantized Hall conductivity correspond to the 
change in density of states when the magnetic flux quanta in the system is increased 
by one and, whereas the quantum number r is the change in density of states when the 
period of the potential is changed so that there is one more unit cell in the system) 15’]. 

For any value of the magnetic flux , the system described by the Hamiltonian ([2]), 
supports only n = 0 solution of Eq. 0 for the quantum numbers a and r. This is 
due to the absence of any gap closing that is essential for topological phase transition 
to states with higher values of cr, r. However, the DE which relates continuously vary¬ 
ing quantities p and (j) with integers cr and r, has some important consequences about 
topological changes in close vicinity of rational values of </>. We now show that the 
infinity of solutions depicted in Eq.Q reside in close proximity to the flux (f> and label 
the fine structure of the butterfly. We illustrate this later in the paper by time-dependent 
perturbation that drives the Hamiltonian ^ periodically. 

We begin by solving the Diophantine equation, using a geometrical approach well 
known in quasicrystal literature - commonly referred as the “Cut and Projection Met hod” \ 11, 
(12). Note that although the explicit solution has been known0, our approach however 
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Figure 2: Geometry associated with the cut and projection method to obtain the solutions of the Diophantine 
equation. In this example, the flux vector is chosen to give </> = 2/3. The set of parallel lines gives the 
solutions for each gap r. Two solutions are shown here, Ti and T 2 . 

illustrates simplicity underlying the number theoretical approach to solve this equation. 
The basic idea is to obtain solutions by going to higher dimensions and the required 
solutions are the projections from two to one dimension. 

We start by defining two vectors: a “flux vector” F and a “topology vector T r as, 

F = (p,q), T r = [ay, r r ] (6) 

The DE is then rewritten as, 


r = F • T r (7) 

This implies that the gap index r is a projection of the topology vector onto the flux 
vector. However, this projection is an integer. This suggests the following scheme to 
obtain the components of T r , namely oy, r r in terms of r as follows. 

As shown in Fig. we consider a two dimensional space with coordinates (x, y). 
A square lattice is defined in this space by considering x and y at integer values. 

A vector F is draw which points to square lattice point (p, q). Then all possible 
solutions of the Diophantine lay in the 2D square lattice, and are contained in a family 
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of parallel lines given by, 


r = (p, q ) • (x, y) (8) 

as indicated in Fig. [2] To find integer solutions, we look at the family of parallel 
lines. These lines are all perpendicular to the vector, 

F ± = (- q,p ) (9) 

which defines the line y = with slope </>,i.e., tana = <p where a is the angle 
between the x axis and F ± . If x is chosen to be an integer, that we associate with a 
Chern number, this will produce a y coordinate, 

y = ( 10 ) 

However, although y is in the family of parallel lines, it does not produce an integer. 
But we find an integer just by taking the floor function of the previous equation, 

[y\ = -\.<t>Pr\ ( 11 ) 

so for each oy, the corresponding r r is given by, 

T r = -L0ovj (12) 

Thus, gaps are labeled by the coordinates of a two dimensional lattice, 

[<T r , T r ] = [oy, |_0(7 r J ] (13) 

Furthermore, by using the identity x = [x\ + {x} to express r r and inserting the 
solution into Eq. ([7]), we obtain that, 


r = q{(/)cr r } (14) 

Notice that care must be taken for negative Chern numbers, since {— x} = 1 — {x} 
for x > 0. We now define the Hull function f as, 

/(0,cr) = {0CT r } (15) 
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which is the filling factor for a Chern number at a given 0. 

This formula can be inverted using the same methodology giving the Chern num¬ 
bers as a function of the gap index, 

° r = I ~ q \^ r + ( 16 > 

The hull function lfl3ll can be viewed as a kind of “skeleton butterfly” plot that en¬ 
codes the topological structure of the Hofstadter spectrum as we explain below. Earlier 
studies have discussed this skeleton in terms of the integrated density of states fT6l . In 
this paper, we use this Hull function along with the numerically obtained butterfly to 
lay out the topological patching of the entire butterfly. Upper and lower graphs in Fig. [3] 
illustrate the relationship between the butterfly graph and its skeleton version obtained 
from the Hull function. We emphasize that although the quantitative analysis of the ac¬ 
tual energies requires a numerical exploration, many features can be obtained using the 
Hull function. As we discuss below, this includes not only the dominant gaps but the 
fine structure associated with them. Fig. [3] shows the filling fractions r/q = /(</>, a) 
as a function of flux 0 for Chern numbers ±1, ±2. Notice that each /(</>, a) is just 
a saw-tooth function cr distinct branches. The intersection of two branches are points 
where distinct Cherns meet. 

Fet us consider two successive gaps in the butterfly landscape, one with gap index 
r = a , Chern cr a and another with gap index r = b , Chern a & emanating from the 
left of the graph and meeting at a certain cj> = p/q. These two gaps will exchange their 
corresponding value of r at the meeting point. From the Diophantine equation with 

(j) 4> + 8(f), we obtain, 

lim_({(<£- t><t>)va} = lim {(<£- S<j))a b }) (17) 

8(f)—^ 0 <50—>-0 

This equation can only be satisfied provided, 

cr a = (Tb + q (18) 

since the fractional part function has period 1 and {x} = x for 0 < x < l,i.e., when 
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applied to this particular case, {(; her} = {{p/q)cr} = ( p/q)cr and the period is q. Thus 
the arguments of the fractional parts in Eq. (17] ) can differ up to multiples of the period. 

This equation, which we will refer as the Chern meeting formula relates the topo¬ 
logical quantum numbers of two swaths of the butterfly that meet at a point. 

We now identify a “central butterfly” and its fine structure as follows (see Fig. <0)- 
The central butterfly is the meeting of the two smallest Cherns (in magnitude ) and the 
value of (j) where they meet forms the center of that butterfly. However, there is a whole 
set of other Cherns that converge at such point. Such Cherns will be identified with the 
fine structure of this central butterfly as we will explain later. 

Having identified the center, which we label as p c /g c , we now ask where are the 
boundaries of this butterfly. A close inspection of the skeleton graph shows that the 
boundaries of the central butterfly are the the closest intersections labeled by the largest 
Chern, larger than the Chern that labels the butterfly. This is illustrated in Fig. [3] 

The precise rules for locating the boundaries of the butterfly specified by a triplet of 
rationals labeling the left boundary, the center and the right boundary , Pc/qc, PR./qR.) 

were discussed in an earlier paper0 which we state as follows. For the central butter¬ 
fly, the Farey neighbors of p c /q c are pl/ qL (left neighbor) and PR/qR( right neighbor) 
where qL < q c and qR < q c . These left and right neighbors of p c /q c define the left and 
the right boundaries of the butterfly. In other words, 


Pc _ PL +PR ( 

q C qL + qR HL UK 

We note that Farey neighbors satisfy the so called “friendly number '’ equation, 



(19) 


( 20 ) 


p c q x - p x q c = ±1 


where x = L,R. 

We now address the question of determining the Chern numbers of every butterfly 
in the butterfly digram, obtaining the topological map of the entire butterfly graph. 
Every butterfly in this graph is characterized by a pair of Chern numbers, which we 
will denote as (cr + , <j_), representing the Chern numbers of the two diagonal wings of 
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Figure 3: The skeleton and the butterfly graphs showing explicitly the Cherns ±1 (red) and ±2 (blue) as 
a function of <f>. Positive Cherns are solid lines and negative Cherns are broken lines. In addition to the 
entire butterfly, we identify three butterflies in the graph, whose centers are shown with black circles and are 
marked in the bottom graph with a trapezoid or a rectangle. The bottom butterfly graph shows these butterflies 
enclosed inside the trapezoids: (2, —2) centered at <j) = 1/4 and (1, —2) centered at 1/3. In addition, there 
is another butterfly enclosed in a rectangle, centered at 1/2. The flux values for three different butterflies are 
also shown with double arrowed lines that are labeled by the filling fraction and the Cherns of the butterfly 
in the upper graph. 
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Figure 4: Three butterflies, all centered at 4> = 3/8 with edges at 1 /3 and 2/5 . The blowup of these 
butterflies (in green) labeled with their Cherns is shown on the right. For off-centered butterflies, thin red 
lines aid the eyes to identify the butterflies. We note that the topological states of the three butterflies are 
linked by cr+ — cr_ = q c = 8. 



Figure 5: Four butterflies, all centered at </> = 2/7 with edges at 1/4 and 1/3. The blowup of these 
butterflies ( in green) labeled with their Cherns is shown on the right. For off-centered butterflies near 
E = 0, thin red lines with black dot as the center aid the eyes to identify the butterflies. We note that the 
topological states of the four butterflies are linked by cr + — <j- = q c = 7. 
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Figure 6: Labeling some of the gaps of the butterfly plot with their Chern numbers. Figure illustrates that 
the higher n-solutions of the Diophantine equation for a given flux reside in close vicinity to that flux. 
Corresponding to certain rational flux values, shown explicitly on the horizontal axis, the color coded dots 
represent gaps whose left and right Cherns are showed inside the parentheses, labeled with matching colors. 

the butterfly. 

The simplest case happens to be the butterflies that are centered on the horizontal 
axis of the butterfly graph where the centers of such butterflies are described by flux 
values Pc/Qc where q c is an even integer. As shown in earlier paper 0, and also fol¬ 
lows from the formula ( [14] ), such butterflies are characterized by a pair Chern numbers 

(a + ,<x_=(f,-f). 

Another simple case of topological structure are the butterflies located at the minimum 
or the maximum values of energy in the butterfly diagram. It is easy to show that these 
butterflies are characterized by Chern pairs (±fe, TQl)- 


In general, as shown in Fig. [4]and Fig. [5] with every Farey triplet {pl/ql^Pc/Qc^Pr/qr) 
that satisfy the Eq. 19 there are q x butterflies where q x = Min(qL , qn). These butter¬ 
flies share the same value of the magnetic flux for their center and the boundaries, but 
are displaced in energy. The Chern number of this family of butterflies are related by 
the solutions of the Chern meeting formula ( Eq. [18]). In other words, these butterflies 
are topologically linked by the equation, 


W+ -cr-l = Qc 


( 21 ) 
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Fig. [4] and Fig. [5]provide examples of possible Chern numbers of such butterflies 
centered at 3/8 and 2/7. Below, we discuss a corollary of the Diophantine equation that 
facilitates the precise rules for determining these solutions and laying out the detailed 
topological map of the butterfly at all values in E and 0. 

To determine the topological map of the hierarchical set of gaps, we study the 
fine structure of the butterfly landscape determining the Chern numbers near a flux 
0o = Po/Qo- To do this, we do a simple “titling” of the flux and p. We substitute in 
the DE 0 = 0o + £0 and p = po + Sp, and the corresponding quantum numbers as 
a = do + Act and r = To + At. Now, taking the limits as £0 and Sp go to zero, we 
obtain, 


4> 0 Aa + At = 0; (22) 

At po 

We will refer this equation as the “corollary ” of the DE equation. 

Since both Act and At are integers and po and qo are relatively prime, the simplest 
solutions of Eq. ( [22] ) are A a = ±nqo and At = where n = 0,1, 2,.... These 

solutions describe the fine structure of the butterfly near a flux 0 O . 

The above corollary determines the entire topological map of the butterfly as de¬ 
scribed in Fig. In addition to the dominant gaps, this plot illustrates the fine 

structure both near the center as well as at the boundaries of every central butterfly in 
the entire butterfly landscape. As seen in the Fig. near 0 = 1/2, we see two cascades 
of gaps , characterized by upper and lower set of Cherns, forming a kind of fountain 
with fountainhead located at E = 0 and E max . 

Having determined the topological structure of the butterfly landscape where the 
gaps have been the focus of our discussion, we now turn our attention to some interest¬ 
ing characteristics of the bands. We show that the fine structure of the gaps discussed 
earlier are rooted in Van Hove singularities that reside at the band center. In the pres¬ 
ence of a magnetic flux p/q, a single band is split into g-bands. The center of each of 
this sub-band exhibits a Van Hove singularity. Therefore, accompanying the hierarchi- 
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Figure 7: (color on line) Density of states (DOS) as a function of the energy for several fluxes. This plot 
illustrates how Van-Hove singularities exist at every band center, irrespective of its location in flux value and 
energy. Furthermore, they form self-similar patterns as for the butterfly centered at 3/8 and with left and 
right boundaries at 1/3 and 2/5. The next butterfly in this hierarchy 0 has its center at 11/30. For odd q, 
the central band is a scaled version of the square lattice, as shown here for 0 = 1/3 and 0 = 1/5 
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Figure 8: (color on line) Illustrating the topological collapse near a Van-Hove singularity. Circles (red) 
show Van Hove singularities where a sequence of cascades of gaps with both positive and negative Cherns 
annihilate. 
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Figure 9: (color on line) By varying the hopping along y-direction periodically, we can introduce a phase 
transition where the gap at p/q = 1/3 and its vicinity is transformed from a Chern-dbl to the Chern =p2 
state. In other words 1 —>• 1 — 3 at 1/3, consistent with the corollary II as ±1 = 1 =p 3. 
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cal set of bands in the butterfly landscape are also a set of Van Hove singularities, seen 
just with an eyeblink, as dark spots - that is, high density points. In fact low magnetic 
flux limit, discussed in the beginning applies to the neighborhood of every rational flux. 

Fig.0 shows the Van Hove singularities at some of band centers in the fragmented 
spectrum in the presence of magnetic field. This clearly implies that Van Hove are in¬ 
tegral part of every band center, irrespective of its location in energy. In other words, 
no matter how fragmented a band is, its center is always a critical point that hosts a Van 
Hove singularity. This is a consequence of the generic existence of saddle points in 
periodic two dimensional systems CD. Consequently, in the case of incommensurate 
flux where bands have zero measure as band width of every band approaches zero, the 
surviving Cantor set or “dust” encodes a fractal set of Van Hove singularities. The fig¬ 
ure shows zoomed versions of the DOS of the butterfly in the interval 1/3 — 2/5. The 
nested set of plots provide an estimate of energy scaling, which turns out to be close to 
10, in agreement with an earlier estimate. ii 


We next address the question of how the topology of the butterfly is influenced by 
the Van Hove singularities. Fig. © shows the topological landscape in the vicinity of 
Van Hove singularities, illustrating what happens to the zero-field Van Hove as the sys¬ 
tem is subjected to a small magnetic flux 0. The magnetic field that fragments the band 
and the resulting cascade of channels or gaps are characterized by positive and negative 
Cherns interlaced as illustrated in the Fig. Again, this type of behavior is present at all 
band centers and Fig. also shows a new Van Hove singularity and the topological 
landscape in its neighborhood at 0 = 1/2. 


Note that the sequence of Cherns near the Van Hove are the higher order solutions 


of the Diophantine equation, described in Eq. ( |22| ). This brings us to an interesting 
point about the importance of these topological states that are crawling around the Van 
Hove. It turns out that under perturbations of the system, these solutions take over, 
replacing the low Chern states with higher Cherns which then become the dominant 
gaps in the system. 
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We illustrate this by perturbing such a systems, to induce quantum phase transi¬ 
tions to topological states with n > 0 given by 0 with dominant gaps character¬ 
ized by higher Chern numbers. We study butterfly spectrum for a periodically kicked 
quantum Hall system fl8l where J y is a periodic function of time t with period-T, 
J y = A S(t/T — n), a system that was recently investigated HU. Readers are re¬ 
ferred to the original paper for various details of the system. 

Figure shows a topological phase transition from Chern ±1 state to Chern ^2 state, 
seen clearly in simple flux values such as 1/3 and 2/3. The new topological state corre¬ 
sponds to the n = 1 solution of the DE as shown in Eq. 0- This can be interpreted as 
a “gap amplification” as the Chem-2 gap after the phase transition has clearly become 
a dominant gap in the butterfly fractal of the driven system. This process where tiny 
gaps of the Hofstadter butterfly can be amplified may provide a possible pathway to 
see fractal aspects of the butterfly and engineer states with large Chern numbers exper¬ 
imentally. 

In summary, the topological characterization of the butterfly landscape is encoded 
in simple rules that determine the topological map of the butterfly at all energy and 
flux scales. The entire hierarchy of miniature butterflies that exists at any scale is ac¬ 
companied by an orderly set of cascades of topological states that attribute topological 
character to the Van Hove singularities. Near these singularities, higher -n solutions of 
the Diophantine equation play a central role. These solutions become dominant under 
perturbations that will induce topological phase transitions to higher Chern states in the 
system. 
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